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Abstract 

In the past few years, a lot of evidences have been accumulated, which indicate that the gluon 
polarization inside the nucleon is likely to be small at least at the low renormalization scales. On the 
other hand, the recent lattice QCD analyses suggest that the net orbital angular momentum carried 
by the quarks is nearly zero. There is also some indication noticed by Brodsky and Gardner based 
on the COMPASS observation of small single-spin asymmetry on the isoscalar deuteron target, that 
the gluon orbital angular momentum inside the nucleon is likely to be small. Naively combining all 
these observations, we are led to a rather embarrassing conclusion that the nucleon constituents 
altogether do not carry enough amount of angular momentum saturating the total nucleon spin. 
We show that this somewhat confused state of affairs can be cleared up only by paying careful 
attention to the scale dependencies of the nucleon spin decomposition. 
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I. INTRODUCTION 



If the intrinsic quark spin carries a little of the total nucleon spin, what carry the rest 
of it ? This is the famous "nucleon spin problem" raised by the EMC measurements nearly 
twenty years ago . In the past few years, there have been several remarkable progresses 
toward the resolution of this long-standing problem. Firstly, a lot of experimental evidences 
have been accumulated, which indicate that the gluonpolarization inside nucleon is likely 
to be small at least at the low renormalization scales Ql-Q]. At the least, it is now widely 
accepted that the f/^(l)-anomaly motivated explanation of the nucleon spin puzzle is disfa- 
vored. Secondly, the quark spin fraction or the net longitudinal quark polarization AE has 
been fairly precisely determined through the high-statistics measurements of deuteron spin 
structure function by the COMPASS [21,0] and HERMES groups j^. According to their 
new analyses, the portion of the nucleon spin coming from the intrinsic quark spin is around 
30%. Putting together these two observations blindly, one might be led to the conclusion 
that the rest of the nucleon spin must be carried by the orbital angular momentum of quarks 
and/or gluons. On the other hand, however, the recent lattice QCD simulations indicate 
that the net orbital angular momentum carried by the quark fields is very small or close to 



zero 



10(1 - 181]. Besides, based on the conjecture on the relation between the Sivers mecha- 



nism and the quark and gluon orbital angular momenta [l9|],[20|], Brodsky and Gardner [21 | 
argued that the small single-spin asymmetry observed by COMPASS collaboration on the 
deuteron target [^j is an indication of small gluon orbital angular momentum inside the 
nucleon. 

Naively combining all the observations above, we might be led to the conclusion that 
the nucleon constituents on aggregate do not carry enough amount of angular momentum 
saturating the total nucleon spin. What's wrong with the above deduction? The purpose 
of the present study is to resolve the apparent paradox above. To clear up this confused 
status of our understanding of the nucleon spin puzzle, we propose to carry out an analysis, 
in which a special care is paid to the fact that the decomposition of the nucleon spin is an 
absolutely scale-dependent idea. What plays a central role in this analysis is Ji's angular 
momentum sum rule, supplemented with some additional knowledge listed below. The first 
is the information obtained from the recent theoretical studies of the isoscalar and isovector 
combinations of the nucleon anomalous gravitomagnetic moments, i?2c^'^(0) and i?2o~'^(0), 
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within the lattice QCD as well as within the chiral quark soliton model (CQSM). The 2nd is 
the empirical information on the momentum fractions carried by the quarks and gluons, as 
well as on the longitudinal quark polarizations. The 3rd is the observation, first made by Ji, 
that the total angular momentum fractions carried by the quarks and gluons obey exactly 
the same evolution equations as the momentum fractions of the quarks and gluons do. 

The plan of the paper is as follows. First, in sect. II, we briefly review main predictions 
of the lattice QCD simulations for generalized form factors and the spin contents of the 
nucleon carried out in the past few years. On the other hand, sect III is devoted to new and 
improved investigation of the corresponding generalized form factors within the framework of 
the chiral quark soliton model (CQSM). Next, in sect IV, armed with the knowledge gained 
in the previous two sections, we try to carry out semi-empirical analysis of the nucleon spin 
contents by paying special attention to their scale dependence. Several concluding remarks 
will then be given in sect.V. 



II. LATTICE QCD PREDICTIONS ON NUCLEON SPIN CONTENTS 

Most theoretical analyses of the nucleon spin contents nowadays heavily relies upon Ji's 



angular momentum sum rule [23|- 26|]. According to it, the total angular momentum carried 



by the quark field with flavor q is given as 

J' = \ [A^o(0) + i?lo(0)] = \ [{xY + BU^)]. (1) 

Here, ^20(0) forward (t — > 0) limit of the generalized Dirac form factor A\Q(t)^ which 

is related to the 2nd moment of the unpolarized spin-non-flip generalized parton distribution 
function (GPD) H'^{x,C,,t). It just reduces to the momentum fraction (x)^ carried by the 
quark with flavor q. On the other hand, -820(0) is the forward limit of the generalized 
Pauli form factor i?|Q(t), which is sometimes called the anomalous gravitomagnetic moment 
(AGM). (More precisely, -B|o(0) the contribution of the quark with flavor q to the nucleon 
AGM.) The quantity -B|o(0) is also related to the 2nd moment of the unpolarized spin- 
flip generalized parton distribution E'^{x,^,t), so that it is in principle measurable through 
the high energy deeply virtual Cornpton scatterings (DVCS) and/or deeply virtual meson 



production (DVMP) processes 23|], 26|. Confining to the two flavor case, for simplicity, we 
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have two independent relations : 



(2) 
(3) 



Since the quark momentum fraction (a;)""^^ and are empirically known fairly well, the 

knowledge of B2q'^{^) and i?2o~°'(0) is essential to extract the total angular momentum J" 
and J'^ carried by the u- and ci-quarks. In fact, these are the quantities of central interest 
in several lattice QCD studies lO] - 18|]. Here, we briefly review the relevant predictions of 
lattice QCD studies on the nucleon spin contents in the past few years. 



We first look into the results on 5^o^^(0) and 5^o"'*(0) reported by the QCDSF Collab 



u—d( 



ration in lOj, HI] some years ago. Their predictions are 



^20 



(0) 



0.102 ±0.113, 



B. 



2V'(0) 



0.566 ±0.113. 



(4) 



(We recall that their simulations were performed in the so-called heavy-pion region with 
rTT-TT ~ (640 — 1070) MeV. The values quoted in (j4]) are those extrapolated to the physical 
pion mass. In practice, however, no strong pion mass dependencies were observed in their 
simulations at this stage.) Combining (jlj) with their predictions on A2q'^{0) and A2q'^{0), 
given by 



A. 



u+d 



0.547 ±0.022, A^^'^iO) 



u—d 



0.253 ±0.022. (5) 



they estimated that 



2J" = 0.74 ±0.12, 2r 



0.08 ±0.08. 



Further combining with their results on the quark polarization, 

Au + Ad = 0.60 ±0.02, Au - Ad = 1.08 ±0.02, 



(6) 



(7) 



they concluded that the net orbital angular momentum (0AM) of the quarks is very small 
or consistent with zero : 

2^n+d ^ 0.06 ±0.14. (8) 
An independent studies of B^^'^^O) and i?2o"'^(0) is reported by the LHPC Collaboration 



121- 



15|: 



^^+^^(0) = -0.09 ±0.03, B^q''{0) = 0.67 ±0.03 



(9) 



Using their previous results for the quark momentum fractions as well as the quark longitu- 



dinal polarizations 



AS = 0.682 ±0.018, (10) 
they also estimated the quark orbital angular momentum to get 

L"" = -0.088 ±0.019, L'^ = 0.036 ± 0.013, (11) 

or 

2iu+d ^ -0.104 ±0.038, 2L"-'^ = - 0.248 ± 0.038, (12) 

Their conclusion at this stage was as follows. Both flavor separately give a rather small 



contribution of the order of 17% (7%) for w-quark (d-quark 



to the nucleon spin, due 
Adding further u and d 



to cancellation in quark momentum fraction, spin and B20 jl5 |. 
contributions give a very small and negative total orbital angular momentum. 

Comparing the results of the two groups, one notices several discrepancies. For instance, 
the central value of the QCDSF prediction for i?2c^'^(0) is small and positive, while the 
corresponding prediction by the LHPC group is small and negative. In spite of these dis- 
crepancies, a main conclusion of the two analyses was common : the net 0AM carried by the 
quarks is very small or consistent with zero. As admitted by themselves, however, a main 
problem of their analyses was that these conclusions were obtained from the simulations 
performed with fairly large pion mass, ranging from 640 MeV to 1070 MeV. 

Very recently, both groups carried out more refined analyses of the nucleon spin contents. 
The simulations were extended to much lower pion mass and the results were further ex- 
trapolated to the physical pion mass with the help of chiral perturbation theory. We first 
overview the main results of the LHPC Collaboration [JLT,]. For the chiral extrapolation, they 
tried several versions of chiral perturbation theory, i.e. covariant baryon chiral perturbation 
theory (BChPT), heavy baryon chiral perturbation theory (HBChPT) with and without the 
A resonance. The results obtained with use of covariant BChPT are 

^^+^^(0) = 0.520 ±0.014, v4^o"'^(0) = 0.157 ±0.006, (13) 

B^+^{0) = -0.094 ±0.050, B^q'^{0) = 0.274 ±0.037, (14) 

which give 

2ju+d ^ 0.426 ±0.052, 2 J" = 0.428 ± 0.032, 2 J'^ = - 0.002 ± 0.032. (15) 



On the other hand, the predictions obtained with the HBChPT without the A resonance 
are given only for the isoscalar quantities : 



A^+'^(0) = 0.485 ±0.014, (16) 

5^+'^(0) = 0.050 ±0.049, (17) 

which give 

2ju+d ^ 0.526 ±0.048. (18) 

One sees that the final answers are fairly sensitive to the ways of chiral extrapolation. In 
particular, i?2o^'^(0), one of our central interest, is slightly negative in the covariant BChPT, 
while it is slightly positive in HBChPT. In either case, combined with their new preliminary 
estimate for the quark spin A'^Q'^(t = 0) = AE"+'^, they reconfirmed their previous conclusion 
that the net quark orbital angular momentum is nearly zero. 



The QCDSF-UKQCD Collaboration also carried out a similar analysis 18|]. Their main 
results are summarized as 

A^^'^iO) = 0.572 ±0.012, A^q'^{0) = 0.198 ±0.008, (19) 
^^+^^(0) = -0.120 ±0.023, S^o"'^(0) = 0.269 ± 0.020. (20) 

We point out that these new results by the QCDSF-UKQCD group changed considerably 
: 'rom the previous QCDSF predictions obtained in the heavy-pion region several years ago 
l^,[ll|. This would mainly be an effect of chiral extrapolation to the physical pion mass. 
Putting aside moderate changes of A^o^'^(O) and A^^^^), the changes of 5^o^'^(0) and B^o'^{0) 
are drastic. First, even the sign is changed for i?2o^'^(0), although the fact, that its absolute 
value is relatively small, is intact. Also drastic is a considerable (more than a factor of two) 
reduction of the magnitude of isovector B^'^{0). (This is also the case for the old and new 



LHPC predictions for i?2o "^(0) Here, we emphasize that this reduction was 



predicted in our theoretical analysis of B2q'^{0) within the chiral quark soliton model [27||. In 
fact, it was shown there that this quantity has a strong pion mass dependence and that the 
lattice QCD predictions obtained in the heavy-pion region has a danger of overestimating it. 
The QCDSF-UKQCD Collaboration also carried out a new estimate of ii(f'^(0) = AS"+'^ 
[isl and obtain 

^j^u+d ^ 0.402 ±0.048. (21) 
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Combining these, they finally obtain an estimate 

2ju+d ^ 0.452 ±0.026, 2L"+'^ = 0.050 ± 0.054. (22) 

Thus, despite some appreciable changes of the predictions for some generalized form factors, 
a common conclusion of the two lattice QCD groups, that the net quark 0AM is small, 
appears to be reconfirmed also by these new analyses. 

Now we have a dilemma. Neither of the intrinsic quark spin, the gluon polarization, nor 
the quark 0AM seems to carry enough amount of angular momentum to saturate the total 
nucleon spin. Does it mean that the rest of the nucleon spin is mostly carried by the gluon 
0AM ? As already mentioned, however, very large gluon 0AM seems to contradict the 
recent claim by Brodsky and Gardner based on the observed small single-spin asymmetry 



on the deuteron target by the COMPASS group [2l|,[22|. In our opinion, this confused 
status arises because we have not paid enough care to the fact that the decomposition of the 
nucleon spin is a highly scale-dependent idea. Later, we shall carry out an analysis, which 
pays more careful attention to the scale dependencies of the nucleon spin decomposition. 

III. CHIRAL QUARK SOLITON MODEL PREDICTIONS 



In a previous paper 27|], we investigated the generalized form factors of the nucleon 



within the framework of the CQSM. A particular emphasis was put there on the pion mass 



dependence of the relevant quantities. (A similar analysis was carried out also in 28|]- 30|. 



See also [3]J, in which the strong pion mass dependence of the net quark polarization AE 
in the chiral region was pointed out.) We discuss here only the predictions on AI^q'^^O), 
S^o^'^(O), A^o^'^(O), and B^q^O), which provide us with enough information for the nucleon 
spin decomposition. 

A largest discrepancy between the predictions of the CQSM and those of the lattice 
QCD simulations was observed for the isovector AGM i?2o~'^(0) of the nucleon [s?], so that 
we will start our discussion with this quantity. Within the framework of the CQSM, or more 
generally in any other low energy models, the forward (t 0) limits of the isovector Pauli 
form factor B'^^'^it) as well as the AGM form factor B^^'^it) are calculated as the difference 
of the standard and generalized Sachs magnetic and electric form factors at t = as (see 
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27l | for more detail) 

5iV'(0) = Git/o^(0) - Gg^(O), (23) 
i?2V'(0) = 4/^(0) - Ggo^(O). (24) 

For completeness, we list below the theoretical expressions for the above quantities within 
the CQSM. The isovector electric form factor in the forward limit, i.e. G^e^I\o) is just 
reduced to the isovector charge of the nucleon, which denotes that 

GgTo'^(O) = 1. (25) 
On the other hand, the isovector gravitomagnetic moment G^'^q^(O) is given as 

^E,20 W - Mm3I [2 ) ^ Em-Er, 

^ ^ m>0,n<0 

X (m||T||n) I — — ("^1 1 1 1'^) + ('^l I 2 ('^ ■ P) 1 1'^) ) (26) 

with Mjv being the nucleon mass. Here, \n) and E^ are the eigenstates and the corresponding 
eigenenergies of the static Dirac Hamiltonian H with the hedgehog mean field, i.e. 

H\n) = En\n), (27) 

where 

H = + pM[cosF{r) + 275 r ■ f sinF(r) ] , (28) 

with M being the dynamical quark mass. The symbols X]n<o Ylm>o stand for the 
summation over all the occupied and unoccupied single-quark eigenstates of H. (The fact 
that G^e11\^) 

is given as a double sum over the single-quark orbitals is connected with 
the fact that it vanishes at the mean-field level and survives only at the first order in the 
collective angular velocity of the soliton.) 

Concerning the isovector magnetic moment G^I[~^q{0) and the corresponding isovector 
gravitomagnetic moment G^^~2q^(0), some comments are in order. In our previous study {2^, 
we have calculated only the leading-order contributions to these quantities and neglected 
the sub leading l/N^ corrections, for simplicity. In the present study, we shall include the 
latters as well. The reason is because a similar 1/Nc correction (or more concretely, the 1st- 
order rotational correction in the collective angular velocity of the soliton) is known to be 
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important for resolving the famous underestimation problem of some isovector observables, 



like the isovector axial-charge, inherent in the hedgehog-type soliton model [32]- 3J]. Taking 
account of this 1st order rotational correction, the isovector magnetic moment of the nucleon 
consists of the leading 0{Q^) term and the subleading 0{Q^) term as 



^M,10 (0) 



G 



(7=1)^0 
MAO 



(0) + G 



(1=1) n'^ 

MAO 



(0), 



(29) 



where 



M, 



N 



G 



(7=1)^1 



MAO 



(0) 



— I 



Nc Y {n\ \ (a; X a) ■ r ||n), (30) 

n<0 

I (mil T \\n) {m\\ (a; x a) x r \\n). (31) 

^ ^ m>0,n<0 ™ 
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Mn f N, 



Similarly, G^go (0) is given as a sum of the 0(f2°) and the 0{Q}) terms : 

<^M,2o (0) = <^M,2o ^ (0) + Gm,^^ (0)' (32) 



where 



G 



{I=1)Q0 
M,20 



G 



(/=i)Qi 

M,20 



(0) 



(0) 



- Nc { {x X ct) ■ T \\n) + (n| | X ■ r | |n) } 

n<0 



(33) 



(m| I T I |n) 



m>0,n<0 
-Em + -En 



Em. Eri 



(m| I (a; X ck) X T I |n) + (m| \ L x t \ \n) > . (34) 



As usual, the above sums over the eigenstates of H can be evaluated with use of the 
discretized momentum basis of Kahana and Ripka 35], 36||. (Some generalization of the 



Kahana-Ripka basis is necessary for the evaluation of the 0{Q^) terms including double 
sums js^].) Now, we are re ady to show the results of our numerical calculation. Similarly 
to the analysis reported in [23], we see the effect of varying the pion mass m^r, by fixing the 
dynamical quark mass M to be 400 MeV. For that purpose, we prepare self-consistent soliton 
solutions for seven values of m^r, i.e. = 0, 100, 200, 300, 400, 500, and 600 MeV, within the 
double-subtraction Pauli-Villars regularization scheme |38|. Favorable physical predictions 
will be obtained by using the value M = 400 MeV and rriTr = 100 MeV, since this set 
gives a self-consistent solution close to the phenomenologically successful one obtained with 
M = 375 MeV and = MeV in the single-subtraction Pauli-Villars regularization scheme 



391]- 4l|. (For the nucleoli mass Mjv, appearing in the above formulas of the generalized 



form factors, the theoretical consistency requires us to use self-consistent soliton masses. 
Otherwise, fundamental conservation laws like the momentum sum rule would be violated. 



See 



27( 1 for the detail. 



TABLE I: The CQSM predictions for the isovector magnetic moment of the nucleon in dependence 
of the pion mass. See the text for more detailed explanation. 



m^(MeV) 


o£r°(o) 


o£r'(o) 







4.12 


1.14 


5.26 


100 


3.41 


1.24 


4.64 


200 


2.89 


1.39 


4.28 


300 


2.69 


1.53 


4.21 


400 


2.67 


1.66 


4.33 


500 


2.72 


1.80 


4.52 


600 


2.73 


1.99 


4.72 



Table [T] shows the theoretical predictions for the isovector magnetic moment of the nu- 
cleon, in dependence of the pion mass m,r- The 2nd and the 3rd columns respectively stand 
for the 0{Q^) and the 0{Q^) contributions, while their sums are shown in the 4th column. 
One can convince that the 1st order rotational correction is very important for this isovector 
observables. With the favorable set of parameters, i.e. m^r = 100 MeV with M = 400 MeV, 
the theory gives jip — ji^ = G'jv^'^Q^(O) — 4.64, which is remarkably close to the empirically 
known isovector magnetic moment of the nucleon : 

ifip-finT'"^ = 4.70589. (35) 

Next, shown in Table [TTl are the predictions of the CQSM for relevant generalized form 
factors in the forward limit, as functions of m^r, which are necessary to evaluate the isovec- 
tor AGM. The 2nd and the 3rd columns of this table respectively stand for the isovector 
gravitomagnetic moment and the gravitoelectric moment, while the 4th column represents 
the leading-order contribution to the isovector AGM of the nucleon. Note that the numbers 
in the 4th column are obtained as the difference of those in the 2nd and the 3rd columns, 
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TABLE II: The CQSM predictions fo the isovector AGM of the nucleon in dependence of the pion 
mass. See the text for more detailed explanation. 



m^(MeV) 




420^(0) 


5^0 ' (0) 


' (0) 


52o"'(0) 


BIq^{^) at Q2 = 4GeV2 





0.361 


0.228 


0.133 


0.272 


0.405 


0.256 


100 


0.392 


0.276 


0.116 


0.342 


0.458 


0.289 


200 


0.452 


0.327 


0.125 


0.429 


0.554 


0.350 


300 


0.519 


0.350 


0.169 


0.491 


0.660 


0.418 


400 


0.579 


0.354 


0.225 


0.534 


0.759 


0.480 


500 


0.640 


0.347 


0.293 


0.567 


0.860 


0.544 


600 


0.716 


0.328 


0.388 


0.600 


0.988 


0.625 



according to the formula, -61/20 (0) = ^M2o (0) ~ ^e2o (0)- These are the predictions 
already given in our previous paper. What is new here is the 5th column, which represent 
the 1st order rotational correction to the isovector AGM of the nucleon. We have already 
seen that the 1st order rotational correction is very important for reproducing the observed 
isovector magnetic moment of the nucleon. Table [TTl shows that the effect of the 1st order 
rotational correction is even more drastic for the isovector AGM of the nucleon. This is 
because the leading-order estimate of the isovector AGM shown in the 3rd column is ob- 
tained as the difference of the two quantities G^M~-^i^) and G'^2o''(0)5 having the same size 
of magnitude, and a sizable cancellation occurs between them. As a consequence, the final 
predictions of the CQSM for the isovector AGM of the nucleon, given in the 7th column are 
nearly a factor of 3 or 4 larger than our previous results neglecting the 1st order rotational 
correction. 

At this stage, one might be interested in a comparison with the predictions of lattice 
QCD. One must be careful here. Different from the anomalous magnetic moment of the 
nucleon, which is scale independent due to the conservation of the electromagnetic current, 
the anomalous gravitomagnetic moment is a scale-dependent quantity. The predictions of 
the lattice QCD simulations corresponds to the renormalization scale of — 4 GeV^, while 
the predictions of the CQSM is thought to correspond to much lower energy scale around 
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= 0.30 GeV^ Fortunately, by making use of Ji's observation that J'^ and {xY obey 
exactly the same evolution equation, we can figure out the scale dependence of B2q'^{Q)- 
(See the next section, for more detail.) From the predictions of the CQSM for B^-'^{Q) 
given in the 6th column, we have estimated the corresponding values at = 4GeV^. The 
results are shown in the 7th column of Table [TTl For the favorable pion mass parameter 
rTT-TT = 100 MeV, our estimate gives B2q'^{'S) ~ 0.289, which should be compared with the 
corresponding prediction _B2o~°'(0) = 0.274 ± 0.037 of the new LHPC lattice simulation JitI . 
and 5^o"'^(0) = 0.269 ±0.020 of the QCDSF-UKQCD one |l8|. One finds that the predictions 



of the CQSM and those of the lattice QCD simulations are now remarkably close to each 
other. This is a welcome result, since it is thought to give a strong support to the reliability 
of the theoretical predictions on the isovector AGM of the nucleon B2q'^{Q). 
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FIG. 1: The pion mass dependence of B^^'^iO) predicted by the CQSM, in comparison with the 
old and new lattice QCD predictions. Both correspond to the scale ~ 4GeV^ 



After obtaining a refined estimate of B2q'^{0) within the framework of the CQSM, we 
revise Fig. 5(b) in our previous paper The filled circles in Fig.l are the CQSM predictions 
of B2q'^(0) corresponding to the scale = 4GeV^, in dependence of the pion mass. The 
corresponding predictions of the QCDSF and LHPC collaborations carried out in the heavy 
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)ion reK ion several years ago are represented by the open squares and the open triangles 



10j,(l5|. On the other hand, the new predictions by the LHPC and QCDSF-UKQCD 



collaborations, extrapolated to the physical pion mass by utilizing the chiral perturbation 



id 



18|. One 



theory, are shown respectively by the open diamond and the filled triangles 
sees that effect of chiral extrapolation is drastic such that the new predictions of the lattice 
QCD are more than a factor of two smaller than the old predictions given in the heavy 
)ion region. This sizable reduction is just consistent with our analysis based on the CQSM 



271]. Now, one can convince that the predictions of the CQSM and the lattice QCD for the 
isovector ACM of the nucleon i?2o~'^(0) is mutually consistent. 



Next, we turn to the discussion of more difficult isoscalar quantities. As shown in [27|,[28|, 
the CQSM predicts that 

^2o^'(0) = {xr"-' = 1, i?2V'(0) = 0. (36) 

It should be noticed that these equalities hold irrespectively of the pion mass within the 
model. The 1st relation is only natural. It simply means that the momentum sum rule is 
saturated by the quark fields alone in this effective quark model, which does not contain 
explicit gluon degrees of freedom. The 2nd relation holds by the similar reason. From Ji's 
angular momentum sum rule, we generally have (in the two flavor case) 

2(J"+'^+J^') = + B^+'^{0) + {xY + BIq{0) = 1. (37) 

If this is combined with the momentum sum rule of QCD, 

+ {xy = 1, (38) 

we are led to a novel identity, 

B^^'iO) + BUO) = 0. (39) 

which dictates that the total nucleon AGM (quark plus gluon contributions) vanishes iden- 
tically. The answer i?2o^'^(0) = is therefore an inevitable conclusion of any effective quark 
model without gluon flelds. In both of the LHPC and QCDSF lattice QCD simulations 
carried out in the heavy^pion region several years ago, the mag nitude of B'^+'^{0) was found 



to be fairly small (ll|,(l2j. Since B^^'^iO) is equal to the difference of 2 7"+^^ and 



the small values of the lattice QCD predictions for i?2o^'^(0) at this point were interpreted to 
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indicate approximate equality of the total angular momentum and linear momentum frac- 



tions of quarks and gluons as advocated by Teryaev several years ago [43|]- 4^. However, 
the recently performed ChPT fits by the LHPC and QCDSF-UKQCD Collaborations ap- 
pears to indicate a sizable bending through the chiral extrapolation in the low pion mass 
region, leading to negative BI^q'^^O) of the order of — 0.1, although one must be very careful 
about the fact that the final conclusion depends on the ways of chiral extrapolation method 

ii,\m\. 

Under such circumstances, it would be fine if we can give some useful constraint on the 
magnitude of B2q'^{0). To this end, we first recall the fact that i?2o^°'(0) is given as the 2nd 
moment of the forward limit of the unpolarized spin-flip GPD E^^'^{x,C,,t) as 

B^+\0) = j a;E"+<^(a;,0,0)cia;. (40) 

It is important to recognize that the 1st moment of the same quantity gives the isoscalar 
magnetic moment of the nucleon up to a factor of 3 : 

5]^+'^(0) = j 0, 0) rfx = = 3(Kf + /t"). (41) 

'D E'^^'^{x, 0, 0) was calculated within the framework of the 



45| . (There is also a calculation for the forward limit of the 



The forward limit of the G 
CQSM by Ossmann et al. 
isovector GPD i5"^'^(x, 0, 0) within the CQSM 46[].) It is given as a sum of the two part, 
i.e. the contribution of A^^^ (= 3) valence quarks and that of the vacuum-polarized Dirac-sea 
quarks as 

i?"+'^(x,0,0) = (x,0,0) + E:+'^(x,0,0). (42) 

An interesting findings there are that the valence quark term turns out to have a similar 
shape as the corresponding valence term f^^i'^ix) of the standard unpolarized PDF, while 
the deformed Dirac-sea contribution has a strong chiral enhancement near x = 0, which is 
antisymmetric with respect to the transformation x ^ — x. (Note that the antisymmetric 
nature of the Dirac-sea contribution to E'^~^'^{x, 0, 0) means that it gives no contribution to 
its 1st moment.) Following the schematic analysis carried out in j47| (see also [^^j), we 
therefore propose to parameterize the characteristic feature of E^~^'^{x, 0, 0) in the following 
simple form : 

E"+'^(x,0,0) = CflYix) - D6'{x), (43) 



14 



with C < 0, and D > 0. With this schematic parameterization, the 1st and the 2nd moment 
sum rules of E'^~^'^[x, 0, 0) become 

j E''+\x,0,0)dx = 3C = 3{kP + k"), (44) 

j xE''+\x,0,0)dx = C J f'iaAx)dx + D. (45) 

Using the observed anomalous magnetic moments of the proton and the neutron, the 1st 
relation gives 

C = {hF + k'^Y^-p = -0.120. (46) 
On the other hand, the 2nd relation gives 



B^iO) = C f:y{x) dx + D. 



(47) 



As a matter of course, in the CQSM, the valence and the vacuum polarization contribu- 
tions in ( H71) exactly cancel each other so that the identity -820^'^ (0) = holds. Such an 
exact cancellation may not happen in real QCD, which contains the gluon fields as well. 
Nonetheless, it is reasonable to expect that the general shape of E^~^'^{x, 0, 0) predicted by 
the CQSM, especially its chiral behavior observed in the small x region, would be preserved 
when going to real QCD, which in turn strongly indicates that there will be no change of 
sign in the contribution of the sea-quark-like component to i?2o^°'(0). We thus conjecture 
that the coefficient D in ( H7j) is at least larger than or equal to 0. Combining this with 
the fact that /;+'^( x) is smaller than 1 (this is because the sea-quark-like component 
also carries some portion of the total momentum fraction of the nucleon), we would then 
conclude from P7|) that the lower limit of i?2o^'^(0) is —0.12 (= (k^ + k^Y^p). In carrying 
out a semi-empirical analysis of the nucleon spin contents in the next section, we therefore 
take a standpoint that the precise value of B2q'^{0) is still uncertain, but it lies most likely 
in the range 

- 0.12 < 5^o^^(0) < 0. (48) 

This is the main theoretical uncertainty in our semi-phenomenological analysis of the nucleon 
spin contents carried out in the next section. 
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IV. SEMI-EMPIRICAL ESTIMATE OF NUCLEON SPIN CONTENTS 



Now, we are ready to start our semi-empirical analysis of the nucleon spin contents. Our 
strategy here is to use empirical information as much as possible, if available. To explain 
our approach more concretely, we start again with Ji's angular momentum sum rule written 
in a slightly more general form : 

1 = JQ + JS, (49) 

with 



jQ 



J' = \ [{xY + 5fo(o)]. 



(50) 



Here, Q denotes the sum of all active quark flavors. {Q = u + d for the two flavor case, and 
Q = u + d + sioi the three flavor case.) To carry out flavor decomposition of the total quark 
angular momentum, we also need another combination of Ji's sum rule : 



J 



u—d 



u — d I TDU — d 

+ ^20 



(0)]. 



(51) 



We emphasize again that the momentum fractions (x)*^, (x)""*^, and (x)^ are all empirically 
well determined. Naturally, these momentum fractions are all scale-dependent quantities. 
A key observation here, first made by Ji 23|,26l], is that J*^ and obey exactly the 
same evolution equations as (x)*^ and (x)^ do. According to him, the underlying reason is 
that forming spatial moment of energy momentum operator does not change short distance 
singularity of the operator. The solution of this (coupled) evolution equation is extremely 
simple at the leading order (LO) : 



2J<3(g2) 
2J^(g') 



3n/ 



16 + 3n/ 
16 



+ 



+ 



2\ 2{16+3nj)/(33-2nj:) 



InQVA^ 

J^g2^^2\ 2{16+3nj)/(33-2n^) p 



2 j^(g^ 

2J\Ql 



3 n 



f 



16 + 3n/ 
16 



16 + 3n/ ' VlnQVAV L ^ 16 + 3n/ 

Particularly interesting here are the asymptotic values in the Q"^ —>■ oo limit 

3n/ „ , 16 



(52) 
(53) 



2 J«(oo) 
Numerically, we obtain 



16 + 3n/' 



2 J^(oo) 



16 + 3n 



/ 



(54) 



2J^(oo) ~ 0.529, 2J^(oo) ~ 0.471, 



(55) 
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for Hf = 6, while 



for rii 



2J^(oo) ~ 0.360, 2J^(oo) ~ 0.640, 



(56) 



In our actual analysis below, we take account of the scale dependencies of the relevant 
quantities by using the known evolution equations at the next-to-leading order (NLO) for 
the momentum fractions, making full use of the fact that and (x)'^ (and also and (x)^) 
obey the same evolution equations. For the sake of completeness, we write down here the 
relevant NLO equations, which we use in the following analysis. The singlet moments 
and (and also (x)'^ and {x)^) evolve as (see, for example, 48|-j5o|) 




A-/2A) 



1 asm-c^siQ' 



2Po 



An 



P RP 



(A+-A-)/2/3o 



P-RP, 



2/?o + A+-A_ 



+ 



(57) 



Here, as{Q'^) is the QCD running coupling constant at the NLO given by 



An 



A lnln(gVA2) 



(58) 



/?oln(QVA2) /?o'ln(QVA2) 
with the choice A = 0.248 GeV, while /3o = ll — |nj and /5i = 102 — with n/ being 
the active number of quark flavor. The quantities R and P± are defined by 



R 



± 



7 



(0) - A= 



(59) 
(60) 



where 7*-'''' and 7'-^^ are the relevant anomalous dimension matrices at the LO and NLO, 



respectively, given by j5l|-|53[| 



7 



(0) 



64/9 -Anf/3 
-64/9 Anf/3 



and 



7 



(1) 



64 
243 



367- 39 'W^f 
(367- 39^/) 



1833 
32 

1833 



(61) 



(62) 
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while X± are the two eigenvalues of the LO anomalous dimension matrix 7^"^ 



On the other hand, the nonsinglet (NS) moments evolve as 48|- SOj 



Jns{Q') 
with 



47r \ 2Po 



0) 



asiQ' 

64 (u 64 



jNs{Qi),m 



iPs = Y> i'nI = ^(367- 39n,). (64) 
Here, it is understood that, for Uf = 3, Jns stands for either of J^^-* = J'^ — J'^ or J^^^ = 

ju^jd_2 js 

Now, we are left with two quantities, i?2o^'^(0) and i?2o '^(0), which are empirically un- 
known yet. Here, one might be tempted to use lattice QCD predictions for those. In our 
opinion, however, blind acceptance of the lattice QCD predictions at the present stage is 
a little dangerous, especially because there seems to be large uncertainties in the process 
of chiral extrapolation. We proceed slightly more cautiously by taking account also the 
information from a phenomenologically successful low energy model of the nucleon, i.e. the 
CQSM. 

After explaining our general strategy, let us now start our semi-phenomenological analysis 
of the nucleon spin contents. We start with the empirical information obtained from the 



MRST2004 as well as the CTEQ5 fits 



5J],(55|. As already emphasized, these two popular 



PDF fits give almost the same quark and gluon momentum fractions below the energy scale 
~ lOGeV^ Ahhough these PFDs are given basically above — 1 GeV^, we try to 
see what happens if we evolve down these fits to lower energy scale as ~ 0.30 GeV^ ~ 
(600 Me V)^, which is understood to be the energy scale of the CQSM. Using the known 



NLO evolution equations for {x)^ and (x)^, together with the MRST2004 predictions 5J], 

(x)'^^ 0.578, (x)^^ 0.422 at Q2=4GeV^ (65) 

we have estimated the scale dependencies of (x)*^ and (x)^ in the range 0.30 GeV^ < < 
4GeV^. The result is shown in Figj2l One sees that the scale dependencies of the quark 
and gluon momentum fractions are fairly strong below ~ IGeV^. At the low energy 
scale around ~ 0.30 GeV^, one finds that the momentum fraction carried by the quarks 
is nearly 80%, while that of the gluons is about 20%. 

As a matter of course, the standard view is that the applicability range of the perturbative 
QCD is at least above 1 GeV, so that one might be a little suspicious of the physical sig- 
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scale dependence of <x> and <x>' 




FIG. 2: The scale dependencies of the quark and gluon momentum fractions, which reproduce the 
MRST fits at = 4GeV^ 

nificance of such "dis-evolution" to low energy scales. Still, we believe it meaningful by the 



ol 



owing reason. Basically, we are following the spirit of PDF fits by Gliick, Reya and Vogt 



48|],[49|. As is well known, these authors start the QCD evolution at the exceptionally low 
energy scales, i.e. Qj^^ ^ 0.23 GeV^ in the leading-order (LO) case, and Q'^^- ~ 0.34 GeV^ in 
the NLO case. They thus found that, even at such low energy scales, they absolutely need 
nonperturbatively (or dynamically) generated sea-quarks, which may be interpreted as the 
effects of meson clouds. We believe such analyses (somewhat nonstandard from the view- 
point of more conservative use of the perturbative QCD) play an important role to connect 
the physics of nonperturbative QCD in the low energy domain and the perturbative QCD 
in the high-energy DIS domain. In fact, we have carried out several theoretical analyses 
based on the GRV spirit. That is, we use the predictions of the CQSM for various PDFs 
as initial-scale distributions given at the low energy scale around 600 MeV. After evolving 
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them with use of the NLO evolution equation, we compare the resultant predictions with the 
corresponding DIS observables with a remarkable success without any other adjustable pa- 
rameters 39|, 40|], 561]- 591]. Then, we shall continue our analysis by accepting the viewpoint 



that the energy scale between 600MeV and 1 GeV is an important region, which connects 
the low energy nonperturbative physics and the high energy perturbative physics of QCD. 
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FIG. 3: The scale dependencies of the quark and gluon angular momentum fractions. The solid 
and dashed curves respectively correspond to the choices i?2o^'^(0) = and i?2o^'^(0) = —0.12. 



Now, we show in Fig|3] our estimate of the scale dependence of the total angular mo- 
mentum fractions carried by the quarks and the gluons. They are obtained in the following 
way. As argued in [23], if the net quark contribution to the nucleon AGM vanishes, i.e. 
B2q'^{0) = 0, we have extremely simple proportionality relations as 



jQ 



(66) 



which was advocated by Teryaev based on the equivalence principle some years ago 42| 



441 ] ■ Very interestingly, these proportionality relations hold independently of the energy 
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scale, since {J'^,J^) and {{x)'^,{x)^) obey the same evolution equations. Thus, the solid 
curves in Figl3]is nothing different from the curves for (x)*^ and (x)^ in FigjJl On the other 
hand, the dashed curves correspond to another extreme, which is obtained by using the value 
B^+'^{0) = -0.12 at Q2 = 0.30 GeVl With this negative value of B^^'^{0), 2 J« becomes 
a little smaller and 2 becomes a little larger as compared with the case -B2o^'^(0) = 0. 
Still, one notices that, at the scale ~ 0.30 GeV^, the quarks carry about 65% of the total 
angular momentum fraction. At the moment, there is a sizable ambiguity in the magnitude 
of i?2c^'^(0), but we believe that the truth lies between the two extreme cases illustrated in 
FigOl (See the discussion at the end of the previous section.) 

Now, the net orbital angular momentum fractions carried by the quarks can be obtained 
by subtracting AS from 2J^. Since the prediction of the CQSM at ^ 5 GeV^ IS re- 
markably close to the central value AS = 0.33 of the recent HERMES analysis, we use 
this HERMES value here. To make the discussion simple, we shall neglect here the scale 



dependence of AS. (In the MS scheme at the NLO, AS is known to have a weak scale 
dependence due to the coupling with Ag. This scale dependence is very weak, however.) 

Shown in Fig H] are 2L^, AS, and 2 as functions of Q^. Here, the solid and dashed 
curves correspond to the cases i?2o^'^(0) = and i?2ci^'^(0) = — 0.12, respectively. First, let us 
look into the case i?2o^'^(0) = 0. In this case, there is a cross over around ~ 0.7 GeV^ ~ 
(840 Me V) 2, where the magnitudes of 2 , AS, and 2 are all approximately equal, 

2L^ ~ AS ~ 2 7^^ ^ 1/3, (67) 

One sees that is a rapidly decreasing function of Q^, so that, as increases beyond this 
cross over energy scale, L*^ becomes less and less important as compared with AS and J^. 
However, the fact that L'^ is a rapidly decreasing function below 1 GeV conversely means 
that it must be very large at the low energy scale around 600 MeV, which we emphasize is 



qualitatively consistent with the picture of the CQSM [37j].|60|. 

Next, we turn to the case i?2o^^(0) = —0.12. In this case, the crossover, where 2 L"^ ~ 
AS ~ 2 ~ 1/3, occurs around the energy scale Q2 ^ 0.30 GeV^ Although the role 
of quark 0AM is less important as compared with the case corresponding to B2q'^{0) = 0, 
it still carries about 1/3 of the total nucleon spin at this low energy scale. We emphasize 
that this is an inevitable conclusion of believing the QCD evolution equation, since it tells 
us that, at this low energy scale, the gluon (spin plus 0AM) carries at most 1/3 of the total 
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FIG. 4: The scale dependencies of the total angular momentum and orbital angular momentum 
of quarks, together with the net longitudinal polarization of quarks. The solid and dashed curves 
respectively correspond to the choices i?2(^'^(0) = and i?20^'^(0) = -0.12 

nucleoli spin, so that what remains to carry the rest 1/3 of the nucleon spin must be the 
quark 0AM. On the other hand, when going to higher energy scale, say at — 4 GeV^, 
corresponding to the renormalization scale of the lattice QCD calculations, one sees that the 
amount of the quark 0AM becomes much smaller. Still, it is seen to carry nearly 20% of 
the total nucleon spin even at — 4 GeV^. One might suspect that this would contradict 
the conclusion of the lattice QCD analyses. Probably, the main cause of discrepancy can 
be traced back to a little overestimation of the net quark polarization AE in the lattice 
QCD. In fact, the results of the QCDSF-UKQCD group for AS is 0.402 ± 0.024 [18], which 
overestimate a little the central value 0.33 of the HERMES analysis 9|], which we have used 
in our semi-empirical analysis here. (One of the reason of a little overestimation of AS in 
the lattice QCD simulations may be attributed to the so-called quenched approximation, 
i.e. the neglect of the disconnected diagrams.) In our opinion, the quark 0AM fraction of 
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the order of 20% is reasonable enough from the following simple consideration. To convince 
it, we recall that the asymptotic value (the value in the ^ oo limit) of the total angular 
momentum fractions of the quarks and the gluons are extracted from the relations (1521) and 
flS5]) . which follows from the fact that {J^ , J^) and {{x)'^,{x)^) obey the same evolution 
equation. With the realistic case of 6 flavors, we have 

2J^(cx)) ~ 0.529, 2J»(oo) ~ 0.471. (68) 

Subtracting AS ~ 0.33, which is thought to be nearly scale- independent, we thus obtain 

2^^3(00) ~ 0.199. (69) 

Since L'^ is a decreasing function of Q^, the magnitude of 2 at the scale ^ 4GeV^ 
must be larger or at least approximately equal to this asymptotic value, which justifies our 
reasoning above. 

So far, we have concentrated on the analysis of the net quark and gluon contribution to 
the nucleon spin and the net quark contribution to the orbital angular momentum. Now, 
we try to make a flavor decomposition of the quark contribution to the nucleon spin and 
orbital angular momentum, which requires the knowledge of the quantity i?2o~'^(0), i.e. the 
isovector nucleon AGM. Since we want to investigate the scale dependencies of the momen- 
tum fractions and the total angular momenta of the quarks and gluons up to Q"^ = 4 GeV^, 
we use again the NLO evolution equation with 3 active flavors, although we assume that 
the strange quarks carry negligible momentum fraction and AGM at the initial low energy 
scale, for simplicity. As initial conditions of evolution, we need the following quantities at 
Q2 = 0.30 GeV^ : 

= = {xy, (70) 

= (x)"-'^, = (71) 

The singlet moments (x)*^ and (x)^ evolve according to the evolution equation flSTl) . Here, 
we use the initial condition 

= 0.785, {x^-'^ = 0.250, (x)^ = 0.215, (72) 

with {x}'^ = 0, since it gives at = 4 GeV^ 

^^y+d+s ^ 57g^ (^^y ^ 421^ (73) 

= 0.552, (x)"-'^ = 0.158, (74) 
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which approximately reproduces the empirical MRST2004 fit at the same scale [5J]. We 
can make a similar analysis also for the total angular momentum of the quarks and gluons, 
because they obey the same evolution equations as the corresponding momentum fractions. 
To proceed, we need initial conditions for the following quantities 

= \ [{xr^' + fi2V'(o)], 
j"-^ = \ [{xr-' + i?2V'(o)], 
J' = I [{^y + Buo)], 

with the general constraint i?2o^'^(0)+i?fQ(0) = 0. Here, we have assumed = at the initial 
scale. For 5^0+^(0), we consider the two cases again, i.e5^o^"'(0) = and S^o^'^(O) = -0.12. 
For i?2o"'^(0), we use the prediction of the CQSM given by 

52^^(0) = 0.458. (75) 

At first sight, the magnitude of the isovector AGM above seems to be fairly larger than the 
corresponding predictions of the lattice QCD given at = 4GeV^. As already mentioned 
in the previous section, however, after taking account of the scale dependence, we find that 
the above CQSM prediction for i?2o~'^(0) is remarkably close to that of the lattice QCD. 

Now that all the necessary conditions are given at the initial scale = 0.30 GeV^, 
let us first try to estimate the total angular momentum fractions at = 4GeV^, which 
corresponds to the renormalization scale of lattice QCD simulations. First, we show the 
results corresponding to the choice i?2o^'^(0) = 0. We have, at = 4GeV^, 

{x)Q = 0.579, = 0.421, (76) 

= 0.552, (x)"-'^ = 0.158, (x)' = 0.028, (77) 



and 



which gives 



2 J« = 0.579, 2 = 0.421, (78) 
2 7"+^^ = 0.552, 2J''-'^ = 0.448, 2 J" = 0.028, (79) 

2J" ~ 0.500, 2 7^^- 0.052. (80) 
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On the other hand, with the choice i?2o^'^(0) = —0.12 at the initial energy scale, we get 

2J^ = 0.519, 2J» = 0.481 (81) 

2ju+d ^ Q_4gg^ 2J"-'^ = 0.448, 2J' = 0.033, (82) 

which gives 

2J" ~ 0.467, 2 7^^- 0.019 (83) 

Depending on the two choices for 5^o^'^(0). i.e. 5^o^'^(0) - -0.12 or B^+'^{0) = 0, we thus 
obtain an estimate, 

2J" ~ 0.46 - 0.50, 2 7^^- 0.02 -0.05, (84) 

which supports the conclusion of the lattice QCD studies that the total angular momentum 
carried by the (i-quarks is nearly zero at least qualitatively. For reference, we show in FigJS] 
the predicted scale dependence of J" and J'^. 

Now, the information on the quark 0AM can be obtained from J", J'^ and J* by subtract- 
ing the corresponding intrinsic spin contributions. Here, we use the empirical information 
provided by the recent HERMES analysis j^, which gives at = 5 GeV^, 

= AS"+'^+" = 0.330 ± 0.011(theor.) ± 0.025(exp.) ± 0.028(evol.), (85) 
gf = AS""'^ = 1.269 ± 0.003, (86) 
gf^ = AS"+'^-2« _ 0.586 ± 0.031, (87) 

Neglecting the error-bars, for simplicity, this gives 

AS" = 0.842, AS"' = -0.427, AS" = -0.085. (88) 

As is well known, due to the conservation of the flavor nonsinglet axial- current, g\ and g\ 
are exactly scale independent. Then, if we neglect very weak scale-dependence of g^^\ all 
of AS", AS*^ and AS'* are thought to be scale independent. Let us first estimate the quark 
0AM at = 4GeV^ Depending on the two choices for 5^o^'^(0), i.e. 5^o^'^(0) = -0.12 
and B^+'^iO) = 0, we obtain 

2L" = -(0.333- 0.300), 2L'^ = 0.489 - 0.522, 2L' = 0.033 -0.028. (89) 

A prominent feature here is that the magnitudes of L" and L'^ are sizably large with the 
opposite sign such that L" < and L'' > 0, which leads to the inequality 
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scale dependence of 2 j" and 2 j' 

1.0 

0.8 

0.6 

0.4 

0.2 

0.0 
-0.2 
-0.4 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

[GeV^] 

FIG. 5: The scale dependencies of J" and J'^. The sohd and dashed curves respectively correspond 
to the choices S^o^'^(O) = and -B^o^'^(O) = - 0.12 

i.e. the isovector dominance of the quark 0AM. (As aheady discussed, the cancellation 
between and L'^ is not so perfect in our semi-phenomenological analysis as compared 
with the lattice QCD predictions.) To understand the physical meaning of the above unique 
feature, we find it instructive to look into the scale dependence of L^~'^ as well as of and 

Shown in Figl6] are scale dependencies of 2L", 2L'^, and 2L"~'^. (Note that the differ- 
ence L""*^ of L" and L'^ does not depend on the choice of -B2o^'^(0).) One clearly sees that 
L""*^ is a decreasing function of Q^. Since L""'^ is negative, this means that IL""*^! is an 
increasing function of Q^. Actually, this somewhat peculiar behavior of can naturally 
be understood from the definitional equation of quark 0AM : 

2L"-'^(g2) = 2J"~'^(g2) _ AS""'^. (91) 
Since J"~'^((5^) is a decreasing function of Q^, while AS"~'^ is Q^-independent, L"~'^((5^) 
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FIG. 6: The scale dependencies of L" and L'^ as well as The solid and dashed curves 

respectively correspond to the choices i?2o^'^(0) = and i?2o^'^(0) = —0.12 

is a decreasing function of . In particular, since 2J"^'^(oo) = 0, as verified from the 
non-singlet (NS) evolution equation (1631) . one finds that the isovector quark 0AM in the 
asymptotic limit oo is solely determined by the isovector axial-charge of the nucleon 

c/5f=^^ = AS"-'^ as 

2L"-'^(oo) = - g^^'^ = -1.269. (92) 

This is really an astonishing observation, since it means that the quark 0AM in the asymp- 
totic limit, at least its isovector combination, is determined solely by the longitudinal 
quark polarization ! Note that, since there is no room for doubt in using the relation 
L"^ = J"^ — ^AS'^ to extract quark 0AM, this mysterious conclusion is an inevitable conse- 
quence of the following two theoretical postulates : 

• the definition of J'^ through Ji's angular momentum sum rule, J'^ = K^)'^ + -^lo ('-')]• 

• the observation that J'^ and (x)'^ obey the same evolution equation. 
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Anyhow, since the net quark 0AM 2L"+'^(Q2) 

is a rapidly decreasing function of Q , we 
can easily understand the feature that is large and negative, while L'^ is large and positive 
above a few GeV scale. It is an interesting open question whether such a large 0AM of the 
u- and (i-quarks with opposite sign can be verified through some direct measurements like 
the single-spin asymmetry of semi-inclusive reactions depending on the Sivers mechanism 
[lol . which is believed to be sensitive to the 0AM of nucleon constituents. 

Now, we attempt to give a complete solution to our first question, i.e. the problem of 
determining the full spin contents of the nucleon. Our answer on the decomposition of 
the nucleon spin into the sum of L*^, | AS, and is already given in FigH] within the 
range 0.30 GeV^ ^ ^ 4.0 GeV^. Still, uncompleted is further decomposition of into 
the sum of Ag and . Unfortunately, this decomposition is not gauge invariant and it 
cannot be done very reliably as compared with the analysis done so far. Still, the following 
qualitative consideration would be of some help to have a rough idea about the complete 
decomposition of the nucleon spin, thereby clarifying the fairly confused situation pointed 
out in Introduction. A basis of the following analysis is the observation that the gluon 
jolarization in the nucleon cannot be very large at least at the low renormalization scales 
7], 8] and HERMES groups j^. As a simplest trial, we therefore assume that the gluon 
polarization A^f is zero, at the low energy model scale around = 0.3 GeV^. Combining 
this with the CQSM prediction AS = 0.35 for the net quark longitudinal polarization, we 
solve the NLO evolution equation for AS and Ag in the standard MS factorization scheme. 

The resultant AS and Ag as functions of are illustrated in FigJT] together with the 
empirical values obtained in the recent NLO analyses by the COMPASS and the HERMES 
groups [ll-j^, as well as the old SMC fit 61|. As repeatedly emphasized, the new COMPASS 



and the HERMES results for AS are remarkably close to the prediction of the CQSM. Also 
noteworthy here is the strong scale dependence of the longitudinal gluon polarization. In 
spite that we have assumed that Ag is zero at the starting energy scale, it grows rapidly 
with increasing Q^. As nicely explained in 621], the growth of the gluon polarization with 
can be traced back to the positive sign of the relevant anomalous dimension S'^q^J^ . The 
positivity of this quantity dictates that the polarized quark is preferred to radiate a gluon 
with helicity parallel to the quark polarization. Since the net quark spin component in 
the proton is clearly positive, it follows that Ag > at least for the gluon perturbatively 
radiated from the quarks. The growth rate of Ag is so fast especially in the relatively low 
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FIG. 7: The scale dependencies of AE and lS.g, obtained as explained in the text, are compared 
with the recent QCD fits by the COMPASS group (open circle and open triangle) and by the 
HERMES group (filled circle). The old SMC result is also shown for reference by the filled square. 

region that its magnitude reaches around (0.3 — 0.4) already at = 3 GeV^, which may 
be compared with the estimate given by the COMPASS group : 



/\g{Q^ = 2,GeV^)coMPAss ^ (0.2-0.3). 



(93) 



It should be emphasized that the gluon polarization of this size is nothing inconsistent with 
the GRSV standard scenario of the polarized PDF fit [63|. (Almost the same viewpoint was 
emphasized also in a recent bag model study of the gluon polarization 6J].) Let us therefore 
proceed further by assuming that our estimate of shown in Fig|7] is not extremely far 
from the reality, which enables us to carry out a decomposition of into /S.g and . 

FigJH] shows the gluon 0AM obtained in the above way, together with 2 and 2 Ag. 
One sees that the gluon 0AM is a rapidly decreasing function of Q^. This feature 
naturally follows since = — Ag and the increasing rate of Ag is much faster than 
that of J^. Very interestingly, the magnitude of in the vicinity — 1 GeV^ turns out 
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FIG. 8: The spin and 0AM decomposition of gluon total angular momentum as a function of Q^. 

to be fairly close to zero. We are not sure whether this can be interpreted as giving a 
support to Brodsky and Gardner's interpretation of the recent COMPASS observation of 
small single-spin asymmetry on the isoscalar deuteron target. 

Anyhow, keeping in mind that the spin decomposition of the nucleon is highly scale 
dependent, our estimate at the scale ~ 4 GeV^ can be summarized as follows. The net 
angular momentum fractions carries by the quarks and the gluons are 2 J'5~ 0.52 — 0.58 
and 2 ~ 0.42 — 0.48. The total angular momentum carried by the quarks can further 
be decomposed into the spin and 0AM parts as AS ~ 0.33 and 2 ~ 0.19 - 0.25. The 
decomposition of into the sum of Ag and \J' is still very ambiguous. But, the standard 
scenario for the evolution of l\g indicates that the gluon total angular momentum of the 
order of 2 J^' ~ 0.42 — 0.48 is a consequence of the cancellation of relatively large and positive 
Agi and negative gluon 0AM with a httle smaller magnitude. 

Finally, we make a short comment on the recent extraction of the quark total angular 
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momentum through the model-dependent GPD analyses of the semi-inclusive reactions. The 
first experimental result for the quark angular momentum was obtained by the HERMES 
Collaboration by studying the hard exclusive vr*^ production on the transversely polarized 
hydrogen target j65i] . Their results, corresponding to the average energy scale ^ 2.5 GeV^, 
is given by 66|,|67(] 

J" + J72.9 = 0.42 ± 0.21 ± 0.06. (94) 

On the other hand, another combination of J" and J'^ was extracted by the JLab Hall A 
Collaboration through the analysis of the DVCS and the Bethe-Heitler processes on the 
neutron and on the deuteron target 68|. Their result, corresponding to the average energy 
scale ~ 1.9 GeV^, is given by 



J'^ + JVS.O = 0.18 ±0.14. 



(95) 



For reference, we show below the corresponding predictions of our semi-phenomenological 
analysis. Depending on the two choices -B2o^'^(0) = — 0.12 and B2q'^{0) = 0, we obtain 



J" + J72.9 = (0.245- 0.268), = 2.5 GeV^ 



(96) 



and 



J'^ + J75.0 = (0.056-0.078), = 1.9 GeV^ (97) 

Clearly, our estimates lie in the allowed ranges of both the HERMES and JLab determina- 
tions of J" and J'^. However, it is also clear that the error-bars of the two determinations 
are still too large to be able to say something definite. 



V. CONCLUDING REMARKS 



After completing our semi-empirical analysis of the nucleon spin contents, we now try to 
answer several questions raised in Introduction. Accepting the observation that the intrinsic 
quark spin carries about 1/3 of the total nucleon spin, what carry the rest of it ? As we 
have shown, the answer depends on the scale of observation in an essential manner. At 
the relatively high energy scale around ~ 4 GeV^, corresponding to the renormalization 
scale of the recent lattice QCD simulations, the quarks and gluons respectively carry about 
(52 — 58) % and (0.42 — 0.48) % of the total nucleon spin. The total angular momentum 
fraction 2 J*^ carried by the quarks can further be decomposed into the spin and 0AM 
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parts as AE ~ 0.33 and 2L'5 ~ 0.19 — 0.25. Our estimate for the quark 0AM appears to 
contradict the conclusion of the lattice QCD studies that the 0AM carried the quarks is 
nearly zero. The cause of discrepancy can mainly be traced back to a little overestimation 
of the net longitudinal quark polarization AE in the lattice QCD simulation. In fact, once 
we accept to use the central value AE = 0.33 given by the recent HERMES fit, the quark 
0AM of the order of 20% at — 4 GeV^ is nothing unreasonable as can be convinced from 
the following two observations made based on the evolution equations of relevant quantities. 
The one is the fact that the asymptotic (Q^ — > oo) value of the net quark 0AM fraction 
is given by 2L^(oo) = 2 J^(oo) - AE ~ 0.529 - 0.33 ~ 0.199 for n/ = 6. The other is 
the fact that L'^ is a decreasing function of . The decomposition of into the sum of 
the spin and the 0AM parts is still very ambiguous. Nonetheless, the standard scenario 
for the evolution of Agf strongly indicates that the total gluon angular momentum of the 
order 2 — 0.42 — 0.48 at = 4 GeV^ is hkely to be a consequence of the cancellation of 
relatively large and positive /S.g and negative gluon 0AM with a little smaller magnitude. 

At the low energy scales of nonperturbative QCD around ~ (0.30 — 0.70) GeV^, we 
get a very different picture on the nucleon spin contents. In these energy scales, the quark 
0AM, the intrinsic quark spin, and the gluon total angular momentum would give roughly 
the same magnitude of contributions to the nucleon spin, i.e. 2 L'^ ~ AE ~ 2 ~ 1/3. 

Also very interesting is the flavor decomposition of the total angular momentum and the 
0AM carried by the quarks. On the basis of Ji's observation that J* and {xY obey the 
same evolution equation, we have shown that the asymptotic limit of the isovector quark 
0AM is solely determined by the isovector axial-charge of the nucleon or the isovector part 
of the longitudinal quark polarization as 2L"~'^(oo) = — g^J~^^ = — AE""'^ = — 1.269, which 
leads to novel isovector dominance of the quark 0AM at the high energy scale. It is an 
interesting open question whether this unique feature of the quark 0AM at high Q'^ can be 
probed through some direct observations in high energy DIS processes. 
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